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Abstract 

The quantization of the transverse motion energy in the continuous potentials 
of atomic strings and planes can take place under passage of fast charged 
particles through crystals. The energy levels for electron moving in axial 
channeling regime in a system of parallel atomic strings (for instance, [110] 
strings of a silicon crystal) are found in this work for the electron energy of 
order of several tens of MeV, when a total number of energy levels becomes 
large (up to several hundreds). High resolution of the spectral method used 
for energy level search has been demonstrated. Hence this method could be 
useful for investigation of quantum chaos problem. 

Keywords: Quantum effects in channeling, Spectral method, Degeneration, 
Tunneling 



1. Introduction 

The motion of a fast charged particle in a crystal near one of crystallo- 
graphic axes or planes is determined mainly by the continuous potential that 
is the potential of a crystal lattice averaged along the axis or plane, near 
which the motion takes place. The longitudinal component of the particle's 
momentum py parallel to the crystallographic axis or plane is conserved in 
such field. So, the problem on the particle's motion in a crystal is reduced 
to the two-dimensional problem of its motion in the transverse plane. The 
finite motion in the potential wells formed by the continuous potentials of 
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atomic axes orplanes is known as axial or planar channeling, respectively 
(see [H, 0, H, 0, m, @, 0] and references therein) . 

The electron motion under axial channeling described by classical equa- 
tion of motion can be both regular and chaotic. A pronounced example of 
chaotic behavior is axial channeling in a continuous potential created by two 
neighboring atomic strings [110] of diamond-like crystal 

On the other side, the quantum effects can manifest themselves during 
channeling. Particularly, the quantization of the transverse motion energy 
can take place (see, e.g., 0, Ch. 7, §53], [s], Sec. II C]). Investigation of the 
chaotic behavior on quantum level needs statistical analysis of large mas- 
sive of energy levels (thousands or more) Q. Many numerical methods for 
searching the transverse motion energy levels as well as other quantum char- 
acteristics of a particle motion in channeling regime had been developed in 
the pioneering papers on quantum approach to channeling phenomenon (a 
good review of them could be found in the book jfjj]). 

The most suitable for quantum chaos studies is the so-called spectral 
method of energy eigenvalues searching [§]; it was successfully used in in- 
vestigation of chaotic phenomena in nuclear physics [8(. In the channeling 



theory the spectral method had been applied for the first time in [10|, [III, [12 
for the case of low electron energies, E ~ 1 MeV, when the number of energy 
levels in the potential well is small. In that series of papers the evolution 
of incident electron wave function during penetration into a crystal was in- 
vestigated, positions and widths of transverse energy bands were found as 
well as momentum distribution of outgoing electrons, and comparison with 
experimental data had been made. 

Note that use of periodic crystal potential in [lli ll~2"| automatically 
leads to formation of band structure of the energy of transverse motion. In 
the paper 13( high resolution of the spectral method was used to trace out 
how the energy levels in the potential well formed by a single atomic plane 
are split due to possibility of tunneling between neighboring planes. For large 
number of such planes this leads to formation of the band structure. 

The main goal of present paper is demonstration of high resolution of the 
spectral method for the axial channeling case in the range of incident electron 
energies of order of several tens of MeV, where the number of energy levels 
becomes large. To simplify the problem, we consider an electron motion in 
the isolated two- well continuous potential formed by two [110] atomic strings 
of silicon crystal, neglecting the crystal periodicity. 
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2. Method 

The spectral method of searching the energy eigenvalues of the quan- 
tum system [i[ is based on the computation of correlation function for the 
time dependent wave functions of the system at the initial and current time 
momenta, ty(x,y,0) and *$?(x, y, t): 

/oo /*oo 
/ **(x,y,0)*(x,y,t)dxdy. (1) 
-oo J —oo 

Fourier transform of this correlation function, 

/oo 
P{t) exp(iEt/h) dt, (2) 
-oo 

contains information about the energy eigenvalues. Indeed, every solution of 
the time-dependent Schrodinger equation 

H*(x,y,t) = ih^r(x,y,t) (3) 
could be expressed as the superposition 

ty(x,y,t) = A y)exp(-iE n t/h) (4) 

of the Hamiltonian's eigenfunctions u n j(x,y), 

Hu n j(x,y) = E n u n j(x,y), 

where the index j is used to distinguish the degenerate states corresponding 
to the energy E n . Computation of the correlation function ([I]) for the wave 
function of the form (j3j) gives 

P(t)= eM-iE n >t/h)A* ntj A n , d , 

/OO POO 
/ Kj( x > y) u w,j'(x, y)dxdy = 
-oo J —oo 

= exp(-iE n d/ h)A* nJ A n/ j^nn^jj' = 



= 2j \A n ,j\ 2 exp(-iE n t/h). (5) 
Fourier transformation of (JSJ) leads to the expression 

P B = 27r/^K/ (6) 



n>3 



We see that the Fourier transformation of the correlation function looks like 
a series of 5-form peaks, positions of which indicate the energy eigenvalues. 

So, the computation of the energy levels for the given system consists of 
the following steps: 

1. Choosing the arbitrary initial wave function ty(x,y,0). The only condi- 
tions of the choice are: 

a) tendency to zero under x, y — > ±oo, necessary for every bound state; 

b) wide spectrum that covers the depth of the potential well; 

c) absence of any symmetry which could lead to the lack of some eigen- 
functions in the superposition PJ (see the discussion in Sec. |3j). 

Asymmetric Gaussian waveform would be a good choice for the most cases. 

2. Numerical integration of the time-dependent Schrodinger equation ([3]) 
with the initial value ^f(x, y, 0) for the discrete series of the time momenta; 
the value of the time step At as well as other computational details are 
discussed in 0, [l3| . 

3. Computation of the integral ([T]) for every discrete time momentum from 
t = to some maximal t = T. Subsequent integration of the obtained 
correlation function P(t) with the exponent in (T5]) is carried out over the 
finite time interval: 

P E = ! P(t)exp(iEt/h)dt. (7) 
Jo 

As a result, we obtain a series of peaks of finite width (inverse proportional 
to T) instead of infinitely narrow 5-like peaks ([6]). 

Note that it is possible to find the energy eigenvalues by simple Fourier 
transformation of the computed wave function fflh at the fixed point (x = 



y = 0, for instance) as it was done in 12]. However, it is a risk to lose some 



eigenvalue if the particular eigenfunction u n j(x,y) is equal to zero in this 
point. 
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3. Results and discussion 



The motion of the fast charged particle in a crystal under small angle ip 
to the crystallographic axis densely packed with atoms could be (with good 
accuracy) described as a motion in the continuous string potential (e.g. the 
potential of the atomic string averaged along its axis) 0,11]. The longitudinal 
(e.g. parallel to the string axis) component of the particle's momentum p\\ 
is conserved in such a field. The motion in the transverse plane will be 
described in this case by the two-dimensional analog of Schrodinger equation 
@, Ch. 7, §53] 

f h 2 1 d 

{ 2E^ 2 + c/(x ' y) J * (x ' y ' t] = ih ^^ V ' *)' (8) 

where V 2 = d 2 /dx 2 + d 2 /dy 2 is the two-dimensional Laplasian operator, and 
the value E^/c 2 (where E\\ = (m 2 c 4 +p|c 2 ) 1//2 ) plays the role of the particle's 
mass. 

The continuous string potential could be approximated by the formula 
@, Ch. 6, §41] 

U l{x ,y) = -U 0l n(l + x2 / y ? +aR2 ), (9) 

where for the [110] string of silicon Uo = 60.0 eV, a = 0.37, (3 = 3.5, 
R = 0.194 A (Thomas- Fermi radius); the least distance between two parallel 
strings is a/4 = 5.431/4 A (where a is the lattice period). So, the continu- 
ous potential, in which the electron's transverse motion takes place, will be 
described by "two-well" function (Fig. [Q 

U(x, y) = Ui(x, y + a/8) + U x {x, y - a/8) (10) 

(neglecting the influence of far-away strings). The finite motion of the elec- 
tron in such potential (corresponding to negative values of the transverse 
motion energy E±) is known as axial channeling (^j. 

To search the transverse motion energy levels in the potential ( flOl) by 
spectral method we have chosen the initial wave function of asymmetric 
Gaussian form: 



V(x,y,0) = — - — { exp 
na x a y 



(x-x ) 2 (y + a/8-yo) 



21 



2a 2 x 2a 2 



+ 
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+ exp 



(x-x ) 2 (y-a/8-y )' 



2*5 2*J 



where a x = 0.05 A, a y = 0.06 A, x = a/35, y = a/45 (high symmetry of 
the initial wave function could lead to the absence of some eigenfunctions in 
the superposition (T4]) and, as a consequence, to the loss of some energy levels, 
see the discussion below and the dotted curve on the lower plot on Fig. [3]). 

Fourier transformations of the correlation functions for different values of 
En are presented in Fig. [2j One can see that at the electron energy increase 
the transverse energy levels shift themselves deeply into the potential well, 
and the total number of levels in the well increases. This is the manifestation 
of general quantum mechanical connection between the particle mass (which 
role is played by the value -Ey/c 2 in our problem, as it was mentioned above) 
and the ground state energy as well as the total number of levels in the 
given potential well. The semiclassical estimation of the number of levels for 
the channeled electron as a function of its En could be found in [3, Ch. 7, 
§53]. Some details of "sucking-up" of new levels into potential well from the 
continuum are investigated in 12 . 

Results of computation of the transverse energy levels for the En = 20 
MeV electron in double fllOp and single potential wells are presented in 
Fig. [3J The logarithm of the absolute value of the Fourier transform of 
the correlation function fl7|) is plotted vs the E± value in the potential well. 
To simplify the comparison, the plots for both single and double wells are 
displayed as mirrors of each other, and the negative value U\(0,a/A) w —4.11 
eV is added to the potential energy (Q to shift the bottom of a single well 
to the level of double one. 

We see that in the range of E± below the saddle point of the double-well 
potential the general arrangement of energy levels is similar to that for a 
single well. However, the splitting of some levels in double-well potential is 
observed. For deep levels the splitting is caused by the breaking of axial 
symmetry for any potential well in the double-well case. As known, the 
quantum states in two-dimensional potential possessing axial symmetry are 
characterized by two quantum numbers: radial n r (that coincides with the 
number of zero points of the radial wave function except ones at the distances 
r = and r — > oo from the center of the field) and projection m of the orbital 
momentum to the field axis of symmetry (see, e.g., the problem 4.7 in [Hj]). 
The states with m = are non-degenerated, and the states with \m\ ^ are 
twice degenerated (positive and negative m correspond to the same energy). 
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Fourier transformation of the correlation functions (J7J) for electron in a single 
axial symmetric potential well is presented by the lower plot of Fig. El 
The solid curve is computed for the initial wave function of a general form, 
and the dotted curve is computed for axial symmetric wave function. Such 
special wave function contains only m = eigenstates, hence the peaks on 
the dotted curve indicate only the positions of non-degenerated energy levels. 
The comparison with the upper plot of Fig. [3] demonstrates that only levels 
with \m\ 7^ are split. 

Another mechanism of splitting is connected with the tunneling between 
two wells. The effect becomes visible for the levels close to the saddle point 
of the potential fllOp . where the potential barrier becomes thin and easy 
penetrative (see Fig. @J. 

4. Conclusion 

The quantum mechanical problem on charged particle motion in oriented 
crystal under the axial channeling regime is considered. The spectral method 
of energy level search is applied to electron channeling in continuous potential 
of single [110] atomic string of silicon crystal as well as for duplet of such 
strings. The method demonstrates good resolution in the case of electron 
energies of order of several tens of MeV when the total number of energy 
levels increases up to several hundreds. 

High resolution is necessary for investigation of statistical properties of 
large massive of energy levels; that properties could be used to study the 
behavior of a quantum system, which classical analog allows the dynamical 
chaos (axial channeling in a two-string system, see [3|, Ch. 6, §43]). The 
spectral method has demonstrated its effectiveness for the similar problem 
in nuclear physics jsj]. 
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y , Angstrom 



Figure 1: Potential energy (fTU|) of an electron in the field of continuous potentials of two 
neighboring atomic strings [110] of a silicon crystal. 
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E , eV 



Figure 2: Fourier transformations of the correlation functions ([7J computed for the set of 
values of E\\ for an electron in the double potential well (fT0|) . Positions of the maxima 
indicate the eigenvalues of transverse motion energy E±. 




Figure 3: Upper plot: Fourier transformation of the correlation function ([7]) for the electron 
of energy E\\ = 20 MeV in the double potential well (fT0|) as a function of transverse motion 
energy E± ; vertical dashed line marks the saddle point level for the potential (fTU)) . Lower 
plot: Fourier transformation of the correlation function for the electron in the single well 
([9]), shifted by the value —4.11 eV, computed for the initial wave function of general form 
(solid curve) and for the axially symmetric initial wave function (dotted curve). 
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